function [poles, zeros]= pairPZ (poles, zeros)
for i=1:length(poles)
    for j=1:length(poles)
        if i~=j
            if abs(real(poles(i))-real(poles(j)))<1 && abs(imag(poles(i))-imag(poles(j)))<1
                poles(j)=conj(poles(i));
            end
        end
    end
end
for i=1:length(zeros)
    for j=1:length(zeros)
        if i~=j
            if abs(real(zeros(i))-real(zeros(j)))<1 && abs(imag(zeros(i))-imag(zeros(j)))<1
                zeros(j)=conj(zeros(i));
            end
        end
    end
end
end